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ABSTRACT 

Numerical simulations have been completed for a 
variety of designs for a 90 deg elbow duct. The 
objective is to identify a design that minimizes the 
dynamic load entering a LOX turbopump located at 
the elbow exit. Designs simulated to date indicate that 
simpler duct geometries result in lower losses. 
Benchmark simulations have verified that the 
compressible flow code used in this study is 
applicable to these incompressible flow simulations. 


NOTATION 

C p - pressure coefficient 
d - pipe diameter 

m - mass flow (lbm/sec) 

P - static pressure (psia) 

P ref - static pressure (psia) one diameter upstream 

of pipe bend 

R - pipe bend radius 

Re - turbulent Reynolds number 

U - streamwise velocity 

Ub - bulk velocity 

0 - bend angle 

p - density 


INTRODUCTION 

Liquid propulsion rocket engines employ ducts to 
route propellant and oxidizer among the various 
turbomachinery components. Duct design, including 
cross sectional shape, area, and amount of bending as 
a function of axial distance, influences the presence 
and extent of secondary flows and flow separation 


near the wall. These phenomena determine the 
pressure losses generated in the duct. Secondary flows 
also influence the velocity distribution and dynamic 
loads experienced by downstream components. The 
present study was performed in order to identify a 
duct geometry that reduces the dynamic loads 
experienced by a LOX turbopump located 
downstream of a 90 deg elbow. 

This paper presents the results of a series of numerical 
simulations for several design geometries of a 90 
degree elbow, a common feature in radial-engine 
turbomachinery. Prior to presenting the LOX duct 
simulations, results are presented for a benchmark 
simulation (water flow through a 90 degree elbow) to 
demonstrate the applicability of compressible flow 
codes to an incompressible flow. Results are then 
presented for several duct geometries, including 
constant as well as varying cross sectional shape and 
area. 

NUMERICAL PROCEDURES 

Ducts do not generally fall into a single class of 
engine components, e.g., ducts are found both in 
turbomachinery and combustion devices. Therefore, 
the three principal codes used to analyze 
turbomachinery and combustion devices at MSFC 
were applied in the current study. This also presented 
the unique opportunity to cross check the codes. 

The governing equations in the computational fluid 
dynamics (CFD) analysis are the time-dependent, 
three-dimensional Reynolds-averaged Navier-Stokes 
equations. Three CFD algorithms were used during 
this study. The first algorithm is Corsair, a density- 
based, time marching, implicit, finite-difference 
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scheme. Corsair was developed for studying 
turbomachinery components. The procedure is third- 
order spatially accurate and second-order temporally 
accurate. The inviscid fluxes are discretized according 
to the scheme developed by Roe [1]. The viscous 
fluxes are calculated using standard central 
differences. An approximate-factorization technique 
is used to compute the time rate changes in the 
primary variables. Newton sub-iterations are used at 
each global time step to increase stability and reduce 
linearization errors. For all cases investigated in this 
study, one Newton sub-iteration was performed at 
each time step. To extend the equations of motion to 
turbulent flows, an eddy viscosity formulation is used. 
The turbulent viscosity is calculated using the two- 
layer Baldwin-Lomax algebraic turbulence model [2]. 
Algebraically generated H-grids are used to discretize 
the remainder of the flow field. The code has been 
parallelized using Message Passing Interface (MPI) 
and OpenMP application program interfaces (API’s) 
to reduce the computation time for large-scale three- 
dimensional simulations. 

The second and third algorithms were developed for 
analyzing combustion devices. The second CFD 
algorithm used for this study is FDNS, a pressure- 
based, finite-difference Navier-Stokes solver. The 
algorithm is second-order accurate and uses the 
extended k-e turbulence model and wall functions. 
The three-dimensional, steady, incompressible 
version was used with a structured grid. 

The third CFD algorithm used for this study is Chem, 
which is a density-based, finite- volume Navier-Stokes 
solver. The algorithm is implemented using the Loci 
framework, which allows implementation issues such 
as parallel processing to be handled transparently to 
the coding of the CFD algorithm [3]. The turbulent 
viscosity is calculated using the decoupled version of 
the baseline two-equation turbulence model [4]. 
Three-dimensional, structured grids were created 
using a highly-automated grid generator. Since Chem 
is strictly a perfect gas code, a method was introduced 
to approximate liquids. A large gauge pressure was 
chosen such that normalized density fluctuations were 
small. The molecular weight was chosen to match the 
density of the liquid at gauge pressure, and the 
transport properties of the liquid were specified. 


BOUNDARY CONDITIONS 

For incompressible inlet flow the mass flow, total 
temperature, and the circumferential and radial flow 
angles are specified as a function of the radius. The 


static pressure is extrapolated from the interior of the 
computational domain. 

For incompressible outflow the circumferential and 
radial flow angles, total pressure, and the total 
temperature are extrapolated from the interior of the 
computational domain. The static pressure is specified 
at mid-span of the computational exit and the pressure 
at all other radial locations at the exit is obtained by 
integrating the equation for radial equilibrium. 
Periodicity is enforced along the outer boundaries of 
the H-grids in the circumferential direction. 

At solid surfaces the velocity is set to zero, the normal 
derivative of the pressure is set to zero, and the 
surfaces are assumed to be adiabatic. 


BENCHMARK TEST CASE 

To confirm the applicability of Corsair and Chem to 
incompressible flows, two benchmark test cases were 
performed for the flow of water through a 90-degree 
elbow [5]. Laser Doppler Velocimetry (LDV) and 
static pressure data were obtained as part of the 
experiments. The experiment was performed for a 
Reynolds number of 43,000 and a mass flow 3.6 
Ibm/sec. 

The Corsair and Chem codes produced similar results, 
therefore the results from the Corsair code only are 
presented for brevity. Two views of the Corsair 
computational grid are shown in Fig. 1. The grid 
contains 385 points in the streamwise direction, 93 
points in the circumferential direction and 41 points in 
the radial direction. Figure 2 shows a comparison of 
the predicted and experimental velocity contours 1.0 
duct diameter downstream of the elbow. The 
predicted and experimental results exhibit good 
agreement. Figure 3 shows a comparison of 
experiment and computed horizontal and vertical 
mean streamwise velocity U/U B in the turbulent flow 
-0.6 diameters upstream of the bend. Prominent 
features include inlet boundary layers of depth ~0.09d 
and large central region of uniform velocity. Figure 4 
shows a comparison of the predicted and experimental 
wall static pressure distributions in the region of the 
elbow. Good agreement exists between the predicted 
and experimental results, suggesting that the Corsair 
compressible flow code is suitable for certain 
incompressible flow simulations. Results for the 
benchmark case have been reproduced using much 
coarser grids, e.g., 94 percent fewer grid points. 
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LOX PUMP INLET DUCT CASES 


Two duct cases were simulated using Corsair and 
FDNS. Figure 5 shows two computational grids 
(Elbow #1 and Elbow #2) used to represent a 90 deg 
elbow at the inlet to a LOX pump. Elbow #1, 
simulated using Corsair, has a constant, circular cross 
section in the axial direction. The inlet of Elbow #2, 
simulated using FDNS, transitions to a triangular 
cross section at midsection of the bend with the apex 
along the outer wall. The triangular cross section 
transitions to a circular cross section with 15 percent 
greater area. Each computational grid contains 
191x93x41 points in the stream wise, circumferential, 
and radial directions, respectively. Sufficient grid 
spacing was chosen near the wall to resolve the 
boundary layer. Figure 6 shows the computed exit 
velocities for Elbow #1 (left) and Elbow #2 (right); 
vectors are aligned in the radial-circumferential plane. 
A secondary flow is clearly evident in both cases, 
similar to the flow distribution first reported by 
Prantdl [6]. The circulation is symmetric in Elbow #1, 
but asymmetric in Elbow #2. Figure 7 shows the 
predicted velocities for Elbow #1 and Elbow #2 (with 
vectors in axial-radial plane). The expansion of the 
flow path near the exit of Elbow #2 causes it to act as 
a diffuser. Hence, the flow encounters an adverse 
pressure gradient, resulting in separation. Pressure 
losses are considerably less for the Elbow #1 design. 


PARAMETRIC STUDIES 


Results from the LOX inlet duct case studies showed 
that duct cross-sectional shape and area influence 
flow conditions at the exit and, therefore, inducer 
performance. A series of parametric studies was 
therefore undertaken to study the respective 
influences of pipe bend radius, cross sectional shape, 
and area. Figure 8 illustrates the geometry parameters. 
The Chem algorithm was used with structured three- 
dimensional grids. Circular, constant-area, straight 
cross-sections of length 1 -diameter and 2-diameter 
were added at the elbow inlet and exit, respectively. 
Water was chosen as the working fluid for 
consistency with the benchmark results. All cases 
were run at a turbulent Reynolds number of 43,000. A 
gauge pressure of 100 psia was chosen which resulted 
in a molecular weight of 3523 kg/kmol. For the 
following results, normalized density fluctuations 
were indeed small, approximately 0.2 percent. Figure 
9 shows an example grid for bend radius, pipe radius, 
and cross sectional shape of 10.0 in, 3 in, and circular, 
respectively. 


Flow conditions in a bent duct may be described as 
the result of a balance between centrifugal and 
pressure gradient forces, expressed as 

A P U 2 m 

d R 

where AP = (P ou ter - Pinner) is the static pressure 
difference between inner and outer walls of the bend. 
Results for the benchmark case show good agreement 
with this relationship (Fig. 4). 

Pressure Loss 


A pressure loss coefficient may be defined as the 
difference in total pressure between inlet and exit of 
the duct, normalized by inlet dynamic pressure: 


dPtotal _ 
tfin 




( 2 ) 


Figure 10 summarizes pressure loss coefficients as a 
function of duct bend radius and duct diameter for 
circular cross sections; the larger the symbol, the 
greater the loss. Several subsets of case studies are 
summarized here. In the first, pipe radius was fixed at 
1.0 in and bend radius was varied as represented by 
the vertical distribution of symbols. Minimal pressure 
losses occur for the shortest bend radii, since 
boundary layer growth (hence blockage) is 
proportional to distance from the inlet. As expected, 
the strength of the secondary flow is inversely 
proportional to the bend radius. In the second subset 
of cases, the ratio between bend radius and pipe 
radius was fixed at 10:1, represented by the diagonal 
distribution of symbols. Losses are relative constant at 
small pipe radius, but decrease significantly for larger 
radius. This behavior is due to two reasons. First, 
blockage by boundary layer growth (for larger bend 
radius) is offset by larger duct diameter. Second, the 
secondary flow is also weaker at larger bend radius, 
hence less blockage. In the third subset, the bend 
radius was fixed at 10.0 in and duct radius was varied, 
represented by the horizontal distribution of symbols. 
For smaller radii, the trend is generally toward 
decreasing pressure loss with increasing duct radius. 


Flow Angle Variation 


To relate the water-based CFD results to those of the 
two LOX duct cases, the pump affinity law was used, 

Q LOX QhIO 

N N ^ 

iy LOX iy H20 
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where Q and N refer to volume flow rate and wheel 
speed, respectively. The absolute flow angle, one pipe 
diameter downstream of the elbow exit, may be 
defined as 


(j> = tan 


\ v n 


(4) 


where V r is the tangential speed of a hypothetical 
pump blade at a given radius, r, and V f is the 
magnitude of the streamwise flow velocity. Figure 1 1 
illustrates absolute flow angle for the case of a bend 
radius and pipe radius of 10 in and 3 in, respectively. 
Figure 12 shows the azimuthal variation in absolute 
flow angle at 10, 50, and 90 percent of pipe radius. 
Maximum angular differences are 0.37, 10.94, and 
11.57 deg, respectively. The corresponding incidence 
angle (i.e., angle at which the flow meets the blade) is 
not constant, except perhaps near the hub of the 
pump. Obviously, this behavior cannot be eliminated 
simply by manufacturing the pump with a different, 
radially-dependent inlet blade angle. 


Figure 13 summarizes maximum difference in flow 
angle at 90 percent span for all of the parametric cases 
shown in Fig. 10; the larger the symbol, the greater 
the angular difference. For constant pipe radius, the 
angular difference is inversely proportional to the 
bend radius. This is consistent with the strength of the 
circulation being proportional to the pressure 
difference between inner and outer bend. For constant 
bend radius, the angular difference is inversely 
proportional to pipe radius. 


DISCUSSION 


Comparison of Figs. 10 and 13 shows that an optimal 
duct design must balance the pressure loss against the 
angular difference, respectively, arising from the 
presence of the secondary flow. Ducts with larger 
bend radius are undesirable owing to large pressure 
loss. Moreover, the largest one or two bend radii 
considered in the present study may not be feasible to 
implement for practical reasons. Ducts with smaller 
pipe radius will generally experience larger flow 
angular difference and larger pressure loss. Therefore, 
the optimal duct design with 90 deg bend should favor 
relatively larger bend radius and larger pipe radius. 


a hypothetical LOX pump. The presence of a 
secondary flow may increase high-cycle fatigue, 
cause flow separation near the pump inlet, cause 
additional flow blockage in the blade passages, and 
lower overall pump efficiency. Large pipe radius is 
preferred to minimize the strength of the secondary 
flow, which leads to variations in flow angle. 
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CONCLUSIONS 

A series of three-dimensional Navier-Stokes 
simulations has been performed to identify an 
optimum design for a 90 deg pipe bend at the inlet of 
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Figure 1. Two views of the computational grid for the benchmark pipe bend case. 
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75 deg of curve 




1-dia downstream 


Figure 2. Comparison of predicted and experimental non-dimensional velocity contours (upper left to 
lower right) at 30, 60, and 75 deg of curve, and 1 duct diameter downstream of the elbow. 
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CFD 


Experiment 


Figure 4. Predicted (left) and experimental (right) wall static pressure distributions in the region of 

the elbow. 
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Elbow Elbow 

#1 #2 

Figure 5. Computational grids for 90 deg elbow cases employing circular (left) and triangular cross 

section. Inlet is from above. 



Elbow #1 Elbow #2 

Figure 6. Comparison of predicted exit velocities for Elbow #1 (left) and Elbow #2 (right); vectors are 

in radial-circumferential plane. 
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Elbow #1 Elbow #2 

Figure 7. Comparison of predicted exit velocities for Elbow #1 (left) and Elbow #2 (right); vectors are 
in axial-radial plane. Expansion of cross section near the exit of Elbow #2 induces flow separation. 
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Figure 9. Two views of a grid used for the pipe parametric study (in this example, bend radius and 

pipe radius of 10 in and 3 in, respectively). 
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0 2 4 6 8 

Pipe Radius (in) 

Figure 10. Map of pressure loss coefficients for ducts with circular cross section. A larger symbol 
corresponds to greater losses, hence a less-desirable feature. 
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Figure 11. Example distribution of absolute flow angle at one diameter downstream from 90 deg pipe 

bend. 



Figure 12. Absolute flow angle as a function of circumferential angle (0 deg, inside) at 10, 50, and 90 

percent of pipe radius, based on Fig. 11. 
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Figure 13. Maximum flow angle variation for CFD cases. A larger symbol corresponds to greater 
angular variation, hence a less-desirable feature. 
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